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Artificially-fabricated  semiconductor  microstructures  have  been  found  to  be  very  im¬ 
portant  for  optical  measurements  and  laser  physics,  and  have  been  the  subject  of  study  for 
the  last  couple  of  decades  because  of  their  possible  device  applications  and  also  because 
of  their  intrinsic  simple  physics.  For  modern  day  technologies,  it  is  very  easy  to  judi¬ 
ciously  tailor  any  particular  device  to  exhibit  a  given  physical  property,  especially  when 
the  carriers  are  confined  in  a  certain  dimension,  such  as  quantum-well  wires  and  dots.  In 
quantum  wells,  the  allowed  intersubband  transition  frequencies  are  generally  on  the  order  of 
infrared  frequencies  for  moderate  widths  and  heights  of  the  wells.  Intersubband  and  intra¬ 
subband  transitions  have  been  extensively  studied  experimently1-3  and  theoretically.4-8 
These  studies  involve  a  variety  of  III-V  semiconductor  heterostructures  such  as  GaAs  and 
AlyGai-yAs,  InGaAS  and  InP,  and  InGaAs  and  InAlAs.  In  most  of  these  cases  the 
barrier  regions  are  doped  randomly  to  populate  the  ground  state,  but  in  some  cases  it  is 
also  possible  to  populate  the  ground  state  optically. 

Theoretically,  a  tailoring  of  intersubband  transition  energies  was  proposed  by  Trzeci- 
akowski  and  McCombe4  by  growing  a  thin  layer  of  AlyGa\-yAs  with  variable  Al  fraction. 
This  extra  barrier  in  the  middle  of  the  well  alters  the  ground-state  energy,  with  the  en¬ 
ergy  of  the  first-excited  state  remaining  essentially  the  same  because  of  the  nature  of  the 
wavefunction  of  this  state  which  vanishes  in  the  middle  of  the  well  for  a  symmetric  quan¬ 
tum  well.  Since  the  wavefunction  of  the  ground  state  is  altered  significantly,  the  oscillator 
strength  decreases  as  the  middle  barrier  height  (Al  fraction  in  AlyGa\-yAs)  increases. 

A  similar  tailoring  is  also  possible  by  doping  the  middle  of  the  well  by  a  positive 
delta-function  potential  (6-potential)5.  The  epitaxial  growth  of  the  6-potential  is  done  by 
the  growth-interruption  technique.9  In  the  process  of  growing  the  epitaxy  layer,  the  donor 
or  acceptor  impurities  could  be  grown  randomly.  These  random  impurities  are  very  often 
the  source  of  the  carriers  in  such  structures.  However,  the  growth  of  the  6-potential  could 
be  achieved  by  depositing  appropriate  impurity  atoms  while  suspending  the  deposition 
of  epitaxy  layers.  By  this  technique  a  finite-width  V  or  A-shaped  6-potentials  of  given 
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strength  can  be  formed.  A  finite- width  5-potential  means  that  V a  is  finite,  with  V  and 
a  individually  also  finite,  where  V  and  a  are  the  depth  or  height  and  width  (i.e.  full 
width  at  half  maximum)  of  the  5-potential,  respectively. 

In  the  present  Letter,  we  propose  a  different  kind  of  structure  for  tailoring  of  the 
intersubband  transitions.  Let  us  suppose  a  quantum  well  consists  of  GaAs  of  width  d  and 
height  Vo  with  walls  extending  to  infinity  and  made  of  AlyGa\-yAs  material,  which  is 
doped  randomly  to  provide  enough  electrons  as  carriers  in  the  ground  state.  A  negative 
5-potential  of  strength,  say  5,  is  doped  at  as  from  both  edges  of  the  well.  In  the  present 
calculation,  we  take  the  wavefunctions  and  their  derivatives  divided  by  the  effective  mass10 
to  be  continuoues  at  the  boundaries.  Although  the  epitaxial  growth  of  the  the  5-potentials 
gives  a  finite  potential,  we  treat  it  as  if  it  is  infinite,  meaning,  ^+(a«)  =  ^~(as)  and 
d'*'  ~  ^ =  +(°« )•  Here  5'  =  2™j~,  where  5  is  the  strength  of  the  5-potential 

in  eVA,  ml  is  the  effective  mass  of  the  electrons  in  the  barrier  regions  in  which  the  negative 
5-potential  is  grown,  and  the  ±  superscripts  to  $  refer  to  the  values  of  the  wavefunctions 
at  the  left  and  right  of  the  5-potential.  The  effective  mass  of  the  electrons  is  calculated 
based  on  the  Al  fraction  in  the  AlyGa\-yAs  layers.11  x  is  the  growth  direction  in  which 
electrons  are  confined. 

The  eigenvalue  condition  for  the  5-doped  single  quantum  well  described  above  is  given 
by 

©(£)  =  2{X2  -  Y2}  cos (kd)  +  {(X2  +  Y2)Z~  +  2XYZ+}  sin(Jtd)  .  (1) 

Here,  X  =  2K  +  6',  Y  =  5'exp(-2 Ka6),  Z*  =  {%&  ±  ££,  K  =  y/2m^rE\  and 
k  =  ■  E  is  the  energy  of  the  electron  and  m*w  is  the  effective  mass  of  the  electrons 

in  the  GaAs  layers.11  Equation  (1)  is  appropriate  for  calculating  the  stationary  states  of 
the  quantum  wells  with  negative  or  positive  5-doped  structures  in  the  barrier  materials  by 
finding  the  zeros  the  of  the  0(E),  whereas  the  virtual  (continuum)  states  of  the  structure 
are  calculated  with  the  help  of  transmission  coefficients  which  are  determined  in  the  usual 


manner.  To  calculate  the  virtual  states,  we  consider  the  electrons  to  be  moving  in  the 
AlyGai-yAs  medium  except  for  the  GaAs  layers  with  a  potential  of  —Vo.  The  amplitude 
of  the  i-th  state  wavefunction  A+  in  the  region  left  of  the  ^-potential  in  the  barrier  left  of 
the  quantum  well  is  given  as 

At'(Ei)  =  -4 K2[6'2  +  26'  cos(ka) {2K  cosh(2K  as)  +  6'  smh(2Kas)} 

+  sm(ka){2K6'Z~  sinh(2 Kas)  +  6'2Z+  cosh(2 Kas)  -  Z~(2K2  -  6'2)}].‘1(2) 

Here  Ei  is  the  eigenvalue  of  the  the  i-th  state  and  K  and  k  are  calculated  for  E  = 
E{.  The  other  coefficients  can  be  calculated  by  matrix  multiplications  which  involve  an 
algebraic  exercise  and  are  not  given  here  for  brevity.  The  wavefunctions  of  all  the  states 
are  normalized  to  unity, 

/  *|*(,)i2  =  l,  (3) 

where  b  =  3(d/2  +  as),  which  is  infinity  in  the  practical  sense  since  the  wavefunctions  are 
essentially  zero  beyond  -b  and  b.  The  dipole  matrix  elements  for  the  transition  from  the 
state  i  to  the  state  j  is  expressed  as 


+-*  j)  =  f  dxtyj(x)  x  9i(x)  . 
J  —  6 


Here  ^*(x)  is  the  complex  conjugate  of  the  wavefunction  for  the  state  j.  We  should  point 
out  that  the  complex  conjugate  of  a  wavefunction  for  the  cases  studied  here  is  immaterial 
since  these  states  are  stationary.  However,  it  may  result  in  complex  dipole  moments  for 
the  transitions  between  stationary  and  virtual  states.  The  oscillator  strength  describing 
the  transition  i  *-*  j  is  defined  as 

f(i  <->  j )  =  2m°-*J  1  z(i  <->  j)|2  ,  (5) 

where  u,j  =  L- ;l  is  the  transition  frequency  between  the  states  i  and  j  and  m0  is  the 
mass  of  the  bare  electron. 


We  considere  a  quantum  well  consisting  of  GaAs  material  of  width  120  A  and  barriers 
consisting  of  AlyGa\-yAs  with  y=0.2237.  Both  barriers  are  doped  with  the  ^-potential  at 
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50  A  from  the  edge  of  the  well.  The  effective  mass  of  the  electrons11  in  the  AlyGa\-yAs  re¬ 
gion  for  y= 0.2237  is  0.08669,  whereas  in  GaAs  it  is  0.067.  The  value  of  y=0.2237  gives  a 
180  meV  potential  height.11  For  the  undoped  structure,  there  Eire  three  states  at  21.7,  85.0 
and  173.3  meV.  Figure  1  shows  the  normalized  wavefunctions  of  the  first  four  stationary 
states  along  with  the  potential  profile  in  eV  for  a  -4  eVA  strength  of  the  6-potential.  The 
situation  of  the  fourth  bound  state  is  discussed  below.  It  is  clear  that  the  probabilities  of 
all  the  states  are  greater  in  the  barrier  regions  in  comparison  with  the  undoped  structure, 
are  mainly  localized  in  the  6-potential  regions,  and  most  of  all,  are  symmetric  around 
the  middle  of  the  barrier  because  the  strength  of  the  6-potentials  in  both  barriers  are  the 
same.  (The  term  symmetric  means  that  the  probabilities  to  find  an  electron  in  both  bar¬ 
rier  regions  are  equal  but  the  parities  of  the  states  are  conserved.)  Another  feature  worth 
noting  at  this  point  is  that  the  third  state  is  mostly  localized  in  fhe  barrier  regions,  i.e., 
it  is  bound  to  the  6-potentials.  The  ratio  of  the  probabilities12  in  the  barrier  regions  for 
the  second  and  first  states  is  enhanced  greatly  for  the  doped  structures,  which  favors  of  an 
effect  called  light-induced  drift  predicted  earlier.6 

As  mentioned  earlier,  there  are  only  three  bound  states  for  the  undoped  structures, 
but  beyond  a  certain  strength  of  the  6-potential,  a  fourth  bound  state  appears.  In  view  of 
the  origin  of  the  fourth  state,  we  have  plotted  the  transmission  coefficients  for  the  doped 
structures  for  the  energies  above  the  surface  of  the  well  for  different  values  of  6  in  Fig.  2(a). 
A  virtual  state  for  the  undoped  structure  exists  at  about  350  meV,  whereas  the  6-potential 
creates  a  virtual  state  at  about  200  meV  and  moves  its  energy  down  as  the  strength  of 
the  negative  6-potential  increases.  This  virtual  state  becomes  stationary  beyond  a  -0.9 
eVA  value  of  the  strength  of  the  6-potential.  All  the  energies  are  measured  from  the 
bottom  of  the  quantum  well. 

The  frequencies  of  the  allowed  transitions  u>ij  described  in  Eq.  (5)  and  calculated  with 
the  help  of  the  Eq.  (1)  are  plotted  in  Fig.  2(b).  u>i2  does  not  change  up  to  about  -3.5 
eVA,  and  then  it  drops  gradually  to  zero  which  means  that  after  a  -3.5  eVA  strength  these 
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two  states  get  closer  and  closer  and  finally  vanish  forever.  wi4  decreases  monotonically 
from  -0.9  eVA  up  to  -3.5  eVA.  For  a  range  of  strength  of  the  6-potential,  does  not 
change  much.  The  range  of  strength  for  which  this  plateau  is  formed  is  important  since 
most  of  the  observable  changes  occur  in  this  range,  and  from  now  on  we  will  refer  this  as 
a  "plateau  range”.  When  the  first  two  states  vanish  gradually,  the  upper  two  states  take 
their  place.  This  interchange  of  these  states  is  not  identical  because  the  6-potential  creates 
a  phase  difference  of  tt  over  the  plateau  range.  Changes  brought  in  to  other  two  frequencies 
u>23  and  u>34  by  the  6-potential  axe  also  shown  in  Fig.  2(b),  which  may  be  important  in 
emission  studies  or  optically  populated  absorption. 

Figure  3  shows  the  dipole  matrix  elements  in  (a)  through  (c)  and  the  oscillator 
strengths  in  (d)  through  (f)  for  the  allowed  transitions  as  a  function  of  the  strength  of 
the  6-potential.  The  matrix  element  z(  1  <->  2),  Fig.  (3a),  remains  constant  at  -30  A  up  to 
the  plateau  range,  increases  in  this  range  and  then  decreases  rapidly.  The  constant  value 
of  z(  1  <->  2)  up  to  -3.5  eVA,  which  is  close  to  the  matrix  element  -21.6  A  for  a  similar 
transition  in  the  infinite  well  of  width  120  A,  indicates  that  changes  brought  in  by  the  6- 
potential  to  the  wavefunctions  is  symmetrical  (ignoring  parity).  The  increase  in  z{\  *->  2) 
is  due  to  the  fact  that  the  second  state  moves  down  whereas  the  first  one  remains  the  same. 
This  range  of  6-potential  strength  is  suitable  for  infrared  detectors.  The  oscillator  strength 
shown  in  Fig.  (3d)  goes  through  a  minimum  in  the  plateau  range,  but  only  the  transitions 
up  to  the  bottom  of  this  tip  are  important,  since  the  frequency  u  12  goes  to  zero  after  this 
point. 

There  is  very  interesting  activity  for  the  transition  1  *-*  4.  2(1  «-»  4)  in  Fig.  (3a)  is 
zero  up  to  the  plateau  range  and  goes  through  a  smooth  maximum  in  this  range.  As  we 
have  discussed  before,  the  fourth  state  has  been  brought  in  from  the  continuum  by  the 
6-potential  but  its  wavefunction  is  essentially  the  same  as  that  of  the  first  one.  It  is  the 
plateau  range  in  which  the  wavefunction  of  the  fourth  state  goes  to  the  wavefunction  of 
the  second  one  with  a  phase  difference  of  tt.  Hence,  in  this  range,  while  the  transition  to 
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the  second  state  can  be  used  to  tailor  the  transition,  a  transition  to  the  fourth  state  is  also 
allowed  and  gives  essentially  a  constant  frequency.  This  transition  to  the  fourth  state  is 
only  allowed  in  the  plateau  range  since  this  state  takes  the  place  of  the  second  one  after 
the  plateau  range.  The  matrix  elements  z( 2  <-+  3),  z( 3  «-►  4)  are  shown  in  Figs.  (3b)  and 
(3c)  and  the  corresponding  oscillator  strengths  in  Figs.  (3e)  and  (3f).  If  either  of  these 
second  and  third  states  are  populated  optically,  a  tailoring  of  the  transition  is  possible  even 
for  a  strength  up  to  the  plateau  range.  Moreover,  these  matrix  elements  z(2  «-►  3)  and 
2(3  «-»  4)  and  corresponding  oscillator  strengths  are  useful  in  emission  studies. 

In  conclusion,  we  have  demonstrated  that  a  tailoring  of  the  optical  transition  for  a 
suitable  frequency  is  possible  by  growing  negative  ^-potentials  in  the  barrier  regions  of  a 
single  quantum  well.  We  have  also  shown  that  beyond  a  certain  value  of  the  strength  of 
the  (5-potential,  the  first  two  states  disappear  and  the  upper  two  take  their  places  with  a 
phase  difference  of  n.  If  there  are  not  two  stationary  state  available  to  take  the  places  of 
the  disappearing  ones,  the  <5-potential  brings  one  down  from  the  continuum. 

This  research  was  supported  by  the  U.S.  Office  of  Naval  Research.  The  calculations 
were  performed  at  the  Pittsburgh  Supercomputing  Center  under  Grant  No.  PHY890020P 
and  at  the  Washington  State  University  Computing  Services  Center. 


8 


References 

1.  B.  F.  Levine,  A.  Y.  Cho,  J.  Walker,  R.  J.  Malik,  D.  A.  Kleinman  and  D.  L.  Sivco, 
Appl.  Phys.  Lett.  52,  1481  (1988). 

2.  M.  Nakayama,  H.  Kuwahara,  H.  Kato  and  K.  Kubota.  Appl.  Phys.  Lett.  51,  1741 
(1987). 

3.  K.  Bajema,  R.  Merlin,  F.-Y.  Juang,  S.-C.  Hong,  J.  Singh  and  P.  K.  Bhattacharya, 
Phys.  Rev.  B  36,  1300  (1987). 

4.  W.  Trzecirkowski  and  B.  D.  McCombe,  Appl.  Phys.  Lett.  55,  891  (1989). 

5.  L.  N.  Pandey,  T.  F.  George  and  M.  L.  Rustgi,  J.  Appl.  Phys.  68,  1933  (1990). 

6.  M.  I.  Stockman,  L.  N.  Pandey  and  T.  F.  George,  Phys.  Rev.  Lett.  65,  3433  (1990). 

7.  M.  I.  Stockman,  L.  S.  Muratov,  L.  N.  Pandey  and  T.  F.  George,  Phys.  Rev.  B  45, 
8550  (1992). 

8.  M.  I.  Stockman,  L.  S.  Muratov,  L.  N.  Pandey  and  T.  F.  George,  Phys.  Lett.  A  163, 
233  (1992). 

9.  E.  F.  Schubert,  Optical  and  Quantum  Electronics,  22,  5241  (1990). 

10.  L.  N.  Pandey,  D.  Sahu  and  T.  F.  George,  Solid  State  Commun.  72,  7  (1989).  For  this 
calculation,  we  choose  b=-l. 

11.  The  barrier  height  and  effective  mass  depend  on  y,  the  fraction  of  Al  in  the  layer 
of  AlGaAs  forming  the  barrier.  The  number  of  monolayers  of  AlGaAs  and  GaAs 
define  the  widths  of  the  barrier  and  well,  respectively.  The  barrier  height  is  calculated 
by  the  formula  V{( eV)  =  0.65(1. 155y  +  0.37 y2)  and  effective  mass  by  =  0.067  + 
0.088y,  where  m0  is  the  mass  of  the  bare  electron.  These  formulae  axe  given  in  detail 
in  a  phenomenological  study  by  H.  J.  Lee,  L.  Y.  Juravel,  J.  C.  Woolley  and  A.  J. 
SpringThorpe,  Phys.  Rev.  B  21,  659  (1980). 


9 


12.  The  total  probabilities  of  the  electrons  in  both  barriers  together  for  each  of  the 
first  four  states  are  0.036,  0.503,  0.98  and  0.611  for  a  -4  eVA  6-potential  strength, 
whereas  these  probabilities  for  the  first  three  states  of  an  undoped  structure  are  0.031, 
0.137,  and  0.556.  Since  the  fourth  state  becomes  stationary  at  a  -0.9  eVA  6-potential 
strength,  the  probability  at  that  strength  for  the  fourth  in  is  0.889. 
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Figure  Captions 

Fig.  1.  Energies  in  eV  and  wavefunctions  of  the  first  four  stationary  states  of  a  quantum  well 
of  width  120  A  and  height  180  meV.  Both  barriers  are  doped  by  negative  ^-potentials 
of  strength  -4  eVA  at  50  A  from  the  edge  of  the  well.  The  potential  profile  and 
stationary  states  are  also  shown  in  eV. 

Fig.  2.  (a):  Transmission  coefficient  els  a  function  of  the  energy  above  the  surface  of  the  well 
for  different  5-potential  strengths. 

(b):  Frequencies  of  the  allowed  transitions  as  a  function  of  the  strength  of  the  5- 
potential. 

Fig.  3.  Matrix  elements  in  (a)  through  (c)  and  corresponding  oscillator  strengths  in  (d) 
through  (f )  of  the  allowed  transitions  as  a  function  of  the  strength  of  the  5-potential. 
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